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Abstract: We discuss the matching between long-distance and short- distance at 
next-to-leading in 1/N C and show how the scheme-dependence from the two- loop 
renormalization group running can be treated. We then use this method to study 
the three 0(p 2 ) terms contributing to non-leptonic kaon decays, namely the usual 
octet and 27-plet derivative terms as well as the weak mass term using the Extended 
Nambu-Jona-Lasinio model as the low energy approximation. We also discuss sub- 
tleties in the momentum routing in the low energy theory and a problem in sepa- 
rating factorizable and non-factorizable contributions from the Qq operator in the 
chiral limit. We update our earlier results on the parameter as well. 
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1. Introduction 

The AI = 1/2 rule in kaon decays has been the subject of very many efforts at 
understanding it, see [1] for a review. We briefly discuss it and a short history of 
attempts to understand it in Section 2. In this paper we attempt to put together 
various approaches that have been done before. The short-distance effects are now 
known to two-loops and the extended Nambu-Jona-Lasinio Model enhanced by us- 
ing Chiral Perturbation Theory whenever possible provides a reasonable basis for 
the long-distance description of hadronic interactions needed. We put the two to- 
gether in a way that treats the scheme dependence correctly. The underlying method, 
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reproducing the results of the short- distance running by an effective theory of ex- 
changes of heavy bosons, which we call X-bosons, is discussed in Section 3.2. The 
low energy model is shortly discussed in Section 3.3. In Section 3.1 we recall the def- 
initions of the off-shell two-point functions that we use here to determine the weak 
non-leptonic couplings. The method here is basically to calculate these two-point 
functions to next-to-leading order in 1/N C , but to all orders in the terms enhanced 
by large logarithms involving M^. We then compare with the Chiral Perturbation 
Theory (CHPT) calculations of the same quantity and in the end we calculate the 
relevant physical matrix elements using CHPT. 

In Section 4 we update our earlier results for -B^[2]. Here we discuss in some 
detail the routing issue in Section 4.1, which is rather non-trivial in the presence of 
neutral X-bosons whose direction is not obvious. This also explains the discrepancies 
of the results for very low /x in the chiral limit of [2] and the results of [3]. We give 
therefore updated numbers and expressions for the main results of [2] here. 

Section 5 contains the same discussion but for the AS* = 1 operators Q\ to Q 6 . 
The current x current operators Qi, Q2, and 1 Q3 are computed at next-to-leading 
(NLO) in 1/N C within the ENJL model. The split in Penguin-like and 5^-like 
contributions is discussed. For Q 5 we cannot simply discuss this split, here the 
correct chiral behaviour is only reproduced after summing both contributions. 

When extending the method to Q 6 one discovers that the factorizable contribu- 
tion from Q e has an infrared divergence in the chiral limit. We discuss this problem 
in Section 5.2 and show how it is cancelled by the non-factorizable contribution. 
This problem might be part of the reason why estimates for the Q 6 operator vary so 
widely. After correcting for this we present also results for the matrix elements of 



Finally we put the numerical results for the long- and short-distances together 
in Section 6 and discuss their stability. We also discuss here the coefficients a, b, and 
c defined earlier by Pich and de Rafael [4]. We recapitulate our main results and 
conclusions in Section 7 



The K — > 7T7T invariant amplitudes can be decomposed into definite isospin quantum 
numbers amplitudes as [A = —iT] 
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2. The AI = 1/2 Rule in K -> tttt 
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Where K s ~ X? + e^ 2 °, tf° (2) = (X° - (+)K°)/v / 2, and CP(X 1 ° (2) ) 
In this paper we are interested in the CP conserving part of K — > 7T7t, so we set 
the small phase in the Standard Model CKM matrix elements and therefore e to 
zero. Above we have included the final state interaction phases 5$ and 5 2 into the 
amplitudes A and A 2 as follows. For the isospin 1/2 amplitude 
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With the measured Ks — > 7r°7r° partial width r o, Kg — > tt + tt partial width 
and K + — > 7r + 7r° partial width r +0 [5], we can calculate the ratio 
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This result is what is called the AI =1/2 rule for kaon decays. 

To understand quantitatively this rule has been one of the permanent issues in 
the literature since the experimental determination. It is by now clear that it is the 
sum of several large contributions both from short distance origin [6, 7] and from 
long distance origin [8, 9, 10] which add constructively to make \Aq\ much larger 
than \A 2 \. 

The lattice QCD community has also spent a large effort on this problem, see 
[11] for some recent reviews. 

Among the long distance enhancements of the |A /A 2 | ratio, the order p 4 chiral 
corrections have been found to be quite important. The CHPT analysis to order p 4 
can be found in [10] and both the counter-terms and the chiral logs to that order can 
be found in [12], the chiral logs were originally calculated in [13]. There are some 
small differences between the two results. The fit of the data to both the order p 4 
K — > 7T7T and K — > nnn counter-terms and chiral logs [10, 14] allowed to extract 2 
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to 0(p 2 ), i.e., around 34 % of the enhancement in the AI 
order p 4 and higher CHPT corrections. 



(2.5) 

1/2 rule is due just to 



2 The fit uncertainties to this result were not quoted in [10, 14]. 
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2.1 CHPT to order p 2 

To order p 2 in CHPT, the amplitudes a and a 2 can be written in terms of two 
couplings, 



8 27 

a o = °o a o 
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a 2 = C G 2 7 ZZ ^r-F Q (m 2 K - ml) , 
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(2.6) 
(2.7) 

S = 5 2 = . (2.8) 

The couplings G 8 and G 27 are two of the 0(p 2 ) AS = 1 couplings. They are defined 
in [12] and can be determined from the 0(p 2 ) amplitudes [14] to be 



with 



and 
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0.48 ±0.06. 



(2.9) 



Here we have only included the error bars from the value of the pion decay constant in 
the chiral limit F = (86 ± 10) MeV, this corresponds to f n = 92.4 MeV. Again there 
are uncertainties from the fit procedure and approximations not quoted in [10, 14]. 
Therefore to 0(p 2 ) 
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To understand the difficulty of the task of reproducing (2.5) it is convenient to make 
an l/N c analysis of the 0(p 2 ) result. At large N c , G 8 = G 27 = 1 and 
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i.e. a factor 11.6 smaller than the QCD result in (2.5) ! Notice that to 0(p 2 ) there 
are no quark mass and therefore no chiral logs corrections to the ratio above. So we 
have to explain one order of magnitude enhancement within QCD in the chiral limit 
with \/N c suppressed corrections. 

Another parametrization which will be useful when studying the AI = 1/2 rule 
is the one introduced by Pich and de Rafael in [4] . In this parametrization 
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27 



a + b, 



G 8 = a + b + -(c-b) 



(2.12) 



The nice feature of this parametrization is that a, b, and c have a one to one corre- 
spondence with the three-different QCD quark-level topologies. The ct-type coupling 
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(a) (b) (c) 



Figure 1: The three types of contributions appearing in the evaluation of matrix-elements 
of operators. Namely, (a) Factorizable, (b) -B^-like, and (c) Penguin-like. 

corresponds to configurations that include the factorizable ones (Figure la). This 
coupling is of order 1 in the large N c limit and has only 1/N% corrections. The 6-type 
coupling corresponds to what we call 5^-like topologies (Figure lb) and is of order 
1/N C . This coupling is related to the value of the parameter in the chiral limit. 

The c-type coupling corresponds to what we call Penguin-like topologies (Figure 
lc) and is also of order 1/N C . So in the large N c limit 

a=l and b = c=0. (2.13) 

The main objective of this paper is the calculation of the 1/N C corrections to (2.11), 
i.e. the couplings b and c. 

The coefficients a, b, and c in [4] were defined in a large N c expansion within 
short-distance QCD, i.e. with quarks and gluons. In the low-energy regime where 
the long-distance part has to be evaluated one however cannot distinguish the 1/N% 
corrections to a from the ones to the coefficients b and c. So for us a takes the large 
iV c value a = 1, b — G27 — 1, and c = (3Gg + 2Cr27)/5 — 1. This definition can be used 
both at long and short-distances and only differs by terms of 0(1/N%) with the one 
in [4]. The definition above has also the advantage that all couplings a, b, and c are 
scale independent. Notice that in the present work the l/iV c 2 are from short-distance 
origin only. 

3. The Technique 

3.1 AS =1 and AS = 2 Two-Point Functions 

The theoretical framework we use to study the strangeness changing transitions in 
one and two units was already introduced in Refs. [2, 12, 15]. The original suggestion 
for this type of method was in [16]. The basic objects are the pseudo-scalar density 
correlators 

IT% 2 ) = i J d 4 xe iq - x {0\T (p it (0)P- ? '(x)e irAS =«) |0> (3.1) 

in the presence of strong interactions. Above, a = 0, 1, 2 stands for \ AS\ = 0, 1, and 
2 transitions and i,j are light quark combinations corresponding to the octet of the 
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lightest pseudo-scalar mesons; 

P w ° (x) = —1= ui'-y^u — di^d , P n+ (x) = di^u ,P K °(x) = [si^d] , 
V2 

P^ + (x) = [s«75m] , P m (x) = -k= ui^u + <& 75 <i — 2si7 5 s . 

V6 L J 

(3.2) 

Here and in the remainder, summation over colour-indices inside brackets is assumed 
unless colour indices have been explicitly indicated. These two-point functions were 
analyzed extensively within CHPT to order p 4 in [12]. In that reference we also 
pointed out how on can obtain information on K — > tttt amplitudes at order p 4 from 
off-shell K — > 7r transitions. 

Now, we want to use the 1/N C technique used in [2, 15] to compute the off-shell 
K — > 7r amplitudes and obtain the relevant counter-terms of order p 2 . See [12], for 
explicit details of which counter-terms of order p 4 we can get and possible ways of 
estimating some couplings we cannot get this way. 

In the large N c limit, there is just one operator in the Standard Model which 
changes strangeness by one-unit 

Q 2 = [s 7 M(i _ k) u ]( x )[u^(1 - l5 ) d ](x) . (3.3) 
After the inclusion of gluonic corrections Q 2 mixes with 

Q 1 = [s 7 "(l - l5 )d](x)[u^(l - 75 )«](a;) (3.4) 
via box-type diagrams (first reference in [6]), and with 

Q 3 ee [*f(l - l5 )d](x) J2 [q^(l - l5 )q}(x) 

q=u,d,s 

Qi = [s^{l-l,)dp]{x) £ lfl^-lM(x) 

q=u,d,s 

Q 5 = - l5 )d](x) £ [gj^ + lB)q]{x) 

q=u,d,s 

Q 6 ee [ S Y(1 - li>)dp\{x) £ [f-y^l + ls)q a ]{x) (3.5) 

q=u,d,s 

via the so-called penguin-type diagrams [6]. Since the numerical importance for the 
issues we want to address here is small and for the sake of simplicity we switch 
off electromagnetic interactions. The operator Q4 is redundant and satisfies Q4 = 
Q2 - Qi + Qz- Under SU(3) L x 811(3);? rotations Q- = Q 2 - Qi, Q3, Qa, Q5, and 
Q & transform as 8 L x l R and only carry AI = 1/2 while Q27 = 3Qi + 2Q 2 — Q3 
transforms as 27 L x 1^ and carries both AI =1/2 and AI = 3/2. 

The Standard Model low energy effective action describing \AS\ = 1 transitions 
can thus be written as 

r A5=1 ee -C A5= iEa(/i) / ^yQi{y) +h.c. (3.6) 
i=i J 
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where C AS=1 = (G F /y/5) V ud V* s . 

There is just one operator changing strangeness by two-units in the Standard 
Model, 

Qas=2 = [s^{l - l5 )d](x)[s^(l - l5 )d](x) (3.7) 

which transforms under SU(3)lxSU(3)# rotations as 27 \ x 1 R . 

The matrix elements of the Qi with % = 1, • • • , 6, and Qas=2 operators depend 
on the renormalization group (RG) scale ji such that physical processes are scale 
independent. 

3.2 The X-Boson Method and Matching 

In this section we explain the basics of how to deal with the resummation of large 
logarithms using the renormalization group and how to do the matching between 
the low energy model and the short- distance evolution inside QCD. The guiding line 
here is the 1/N C expansion. 

Let us first explain the philosophy in the case of photon non-leptonic processes 
[15, 17, 18]. The basic electromagnetic (EM) non-leptonic interaction is given by 

CEM = ^f / / d4 ^ / d 4 ye^- y) ^-r Had (x) r Had (y) . (3.8) 

Here we used the Feynman gauge, for a discussion of the gauge dependence see [15], 
= (qQ^q), q T = (u,d,s) and Q is a 3 x 3 diagonal matrix collecting the light 
quark electric charges. The integral over r 2 we rotate into Euclidean space and split 
into a long and a short distance piece, 

J d A r E = J dQ (jT d\r E \ \r E \ 3 + J™ d\r E \ \r E \ 3 ^j . (3.9) 

The long distance piece we evaluate in an appropriate low-energy model, CHPT[18], 
ENJL[15] or using other hadronic models [17]. The short- distance part can be eval- 
uated using the operator product expansion (OPE) and the matrix-elements of the 
resulting operators can be evaluated to the leading non-trivial order in 1/N C using 
the same hadronic low-energy hadronic model as for the long-distance part. 

This procedure works extremely well in the case of internal photon exchange. 
The problem is that in weak decays there are large logarithms present of the type 
ln(Mw I ' Hl) I 'N c which make the 1/N C expansion of questionable validity. The solution 
to this problem at one-loop order was presented in [2] where we showed that the 
integral in (3.9) satisfied the same equation as the one-loop evolution equation. This 
method was very nice for Bk and can also be applied to the AS = 1 transitions. 

Here we will give an alternative description of the method used there that will be 
extendable in a relatively straightforward way to the two-loop renormalization group 
calculations. The precise definition and calculations we defer to a future calculation. 
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We start at the scale M\y where we replace the exchange of W and top quark in 
the full theory with higher dimensional operators using the OPE in an effective theory 
where these heavy particles have been integrated out. So at a scale //# M w we 
need the matching conditions between the full theory and the effective one. As usual 
we get them by setting the matrix elements between external states of light particles, 
i.e. the remaining quarks and gluons, in transition amplitudes with W boson and 
top quark exchanges equal to those of the relevant operators in the effective theory. 

Step 1: at pi H w M w : (2\(W, top-exchange) Full \l) = (2| ]T C t (fi H )Q t \l) . (3.10) 

i 

We then proceed by using the renormalization group to run down from ji H to 
Hi, below the charm quark mass where we have an effective theory with gluons and 
the three lightest quark flavours. At each heavy particle threshold crossed new 
matching conditions between the two effective field theories (with and without the 
heavy particles being integrated out) have to be set, this is done completely within 
perturbative QCD, see e.g. [19]. So that 

Step 2: from pt H to pt L (2| ]T C t (fi H ) Q t \l) — > (2| £ C> L ) Qj\l) . (3.11) 

i j 

At Step 3 we again introduce a new effective field theory which reproduces the 
physics of the operators Qj below fi L by the exchange of heavy Xj-bosons with 
couplings Qi. Again we need to set matching conditions 

Step 3: at \i L : (2|(X,,— exchange) |1) = (2\J2C j (fi L )Q j \l) . (3.12) 

j 

Here the matching means that the left hand side should be evaluated in an operator 
product expansion in M Xl 

The right hand side matrix elements in (3.12) can be evaluated completely within 
perturbative QCD and therefore all the dependence on the renormalization scheme 
and the choice of the basis Qj and of evanescent operators disappears in this step. 
This procedures fixes the gi couplings as functions of the chosen masses and the 
matrix elements {2\J2j Cj(fi L )Qj\l) which are scheme independent. Depending on 
the order to which we decide to calculate in the effective theory, g { will depend on 
additional terms that can be fully determined within the effective theory with heavy 
Xi bosons. 

As an example, let us use the effective field theory with two-loop accuracy for 
the running between scales \xu and fi^ and calculations at next-to-leading order in 
1/N C within the heavy boson effective theory. The term Ci(/j,l) Qi is reproduced 
in the Xi effective field theory by the exchange of a heavy enough vector-boson X± 
with couplings 

AT {gi [s 7 „(l - 75 )d] + g[ [u^(l - l5 )u}} + h.c. (3.13) 
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(a) 



(b) 



Figure 2: The reproduction of the operator Q\ by the exchange of a neutral boson X\. 

The X\ boson has only AS* = 1 components. This is shown pictorially in Fig. 2. The 
scale fx l should be high enough to use perturbation theory. We have the following 
matching conditions (3.12) in this case (we assume that Q\ only has multiplicative 
renormalization for simplicity) 

) =ClW (l + ^r,). (3.14) 

The r\ term cancels the scheme dependence of the two- loop Wilson coefficient Ci(/xl). 
Notice that we can choose independently any regularization scheme on the left and 
right hand sides. In the present work we will use the NDR (naive dimensional 
regularization) two- loop running between fx w and lil- All the large logarithms of 
the type \n(M w / are absorbed in the couplings of the Xj boson in a scheme 
independent way. 

Now we come to Step 4. Assume we want to calculate K° — > tt° matrix element 
in the Standard Model. Since we have included the effect of all the large logarithms 
between Mw and /xl in the gi couplings, we can now apply the same procedure 
explained at the beginning of this section for the photon exchange case [15, 17, 18] 
and remain at next-to-leading order in 1/N C . This we do now for the effective three- 
flavour field theory with heavy massive Xi bosons. So we split the integral over \te\ 
into a long distance piece (between and fx) and a short distance piece (between fx and 
oo) as in (3.9). When evaluating the second term in (3.9) we will find precisely the 
correct logarithmic dependence on M Xl to cancel the one in (3.14). The presentation 
of the scheme dependent constants r\ and f\ for AS = 1 and AS* = 2 is deferred to 
a future publication. 

We then require some matching window in fx along the lines explained in [2] 
between these two pieces. We will use the framework described above to calculate 
AS = 1 and AS = 2 two-point functions and defer the full discussion about this 
procedure to a future publication. In practice we will also choose fx = fx^. 

The same procedure can in principle be used in lattice gauge theory calculations 
where one can then include the X r bosons explicitly in the lattice regularized theory 
or equivalently work with the corresponding non-local operators. 



9i 9i (-, atsifJ-L) 



di In 



M Xl 
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3.3 The Low-Energy Model 

The low-energy model we use here is the extended Nambu-Jona-Lasinio model. It 
consists out of the free lagrangian for the quarks with point-like four-quark cou- 
plings added. This model has the correct chiral structure and spontaneously breaks 
chiral symmetry. It includes a surprisingly large amount of the observed low en- 
ergy hadronic phenomenology. We refer to the review articles [20] and the previous 
papers where we have discussed the various aspects of the ENJL model used here 
[2, 21, 22, 23, 24]. A short overview of the advantages and disadvantages can be 
found in [15] Section 3.2.1. 

It is well known however that it doesn't confine and doesn't have the correct 
momenta dependence at large N c in some cases. These two issues were treated in 
[25] were a low energy model correcting the wrong momenta dependence at large N c 
was presented. 

The bad high energy behaviour of ENJL two-point functions produces some 
unphysical cut-off dependence. In this work we try to smear out this bad behaviour 
as follows. For the fitting procedure we only use points with small values of all 
momenta and always Euclidean. We also keep only the few first terms in the fit to a 
polynomial ( of order six at most) which are therefore not extremely sensitive to the 
bad high energy behaviour of the ENJL model. The model in [25] gives very good 
perspectives that this unphysical behaviour can be eliminated to a large extent, see 
for instance the recent work in [26], and would provide a natural extension of this 
work. 

4. AS = 2 Transitions: Long Distance 

In this section we apply the technique to AS = 2 transitions. These transitions were 
already studied in [2] using the same model for the low energy contributions, there 
are however differences in the routing of the momenta with respect to the one we 
took in [2]. See the next section for a discussion of this issue. 

We study the two-point function U K K °(q 2 ) in the presence of strong interactions 
as defined in (3.1). The operators in La,s=2 are replaced by an X boson coupling to 
[s7 M (l — 75)d](x) currents as described in Section 3.2. 

We evaluate the two-point function then as a function of ji for various values 
of q 2 and masses and this allows us to extract the relevant couplings in CHPT. We 
restrict ourselves here to the 0(p 2 ) coefficient G 2 7 and the actual value of B K . 

4.1 The Routing Issue 

In this section we would like to explain why our present results on B K differ from 
those presented in [2] even though we use the same method and the same model. 
At the same time this will explain the difference between the result from Section 4 
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(a) 



^g) ® «wvw 

(b) 
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■0 
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(d) 







<8> 




<8> 



(c) 



(f) 



Figure 3: Chiral Perturbation Theory contributions to HAS=a(q 2 )- (a) Lowest order. 
(b)-(f) Higher order non-factorizable. The full lines are mesons. The zig-zag line is the 
X-boson. 

in [2] for G27 and the one from [3]. Both papers use the method of [17] and [27] to 
identify the cut-off scale used to identify with the short- distance evolution and we 
have several times checked the calculations in both papers and found no errors in 
either. We will present the discussion here in the case where the low energy model 
used is CHPT to simplify the discussion. 

The source of the difference turned out to be more subtle. In [2] the choice of 
momentum for the X-boson was made to be r + q where q is the momentum going 
through the two-point function defined in (3.1) and r is the loop integration variable. 
This particular choice was done in order to have the lowest order always non-zero, 
even if the range of momenta in r integrated over was such that |r 2 | < \q 2 \. We had 
also always chosen the direction of r + q through the X boson such that the internal 
propagator appearing in diagram (b) of Fig. 3 had momentum r. Since the X in 
that case was a neutral gauge boson this was a natural choice. It turns out however 
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(rX- (rX- (r+2q) 




Figure 4: The routing for the AS = 2 operator enforced by CPS symmetry, (a) Routing 
used in [2], (b)+(c) The correct routing as it should have been used. 



that in the presence of a cut-off some of the contributions obtained with this routing 
do not have the correct CPS symmetry. This symmetry imposes that some of the 
contributions have to have the internal propagator in Fig. 3 with momentum r + 2q 
instead of r. The precise change has been depicted in Fig. 4. The momentum flow as 
depicted in (a) should be replaced by the sum of (b) and (c). This doesn't affect the 
coefficients of the chiral logarithms. Therefore one can use any routing when using 
regularization which doesn't see analytic dependence on the cut-off. Unfortunately, 
this bad routing was actually causing most of the bad behaviour for Bx(n) for high 
values of /i in Table 1 of [2] and the difference with the result for G 2 7 of [2] and [3] . In 
fact, using the background field method as in [3] the CPS symmetry is automatically 
satisfied at order p 2 with any routing. 

We have now corrected for this problem and obtain a much more reasonable 
matching between long-distance contributions and the short-distance contributions. 
Nevertheless, it turns out that the range of values chosen for fi in [2] to make the 
predictions was not very much affected by the routing problem explained above. The 
results we now obtain are much more stable numerically and in the same ranges as 
the ones quoted in [2]. We also agree with the result in [3] for G 2 7(/i) obtained from 
lowest order CHPT, 

g ^ = 1 Sfz- < 41 > 

Here and in what follows, the ji dependent G$(fi), G' 8 (fi), G27(/i), and Bk(^) cou- 
plings stand for the long-distance contributions to those couplings, i.e. with [1 + 
MfJ^MrijlCjifj,) = I. 

4.2 CHPT Results 

Here we update Section 4 of [2] to correct for the routing problem. The non- 
factorizable contribution to Yi K K °(q 2 ) is given by the diagrams in Figure 3 and 
is: 

-8B 2 F 2 j rn d 4 r E r\q 2 E f» d 4 r E r\ 




(2n)Hr 2 E + m 2 K ) 2 J (2rr)4 r\ + m\ 
1 



(r E + 2q E ) 



+ 



{rE + qE? + ml (r E + q E ) 2 + 2m 2 K - ml 
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- 


- 


- 


- 


0.4 


0.737 


0.552 


0.776 


- 


- 


- 


- 


0.5 


0.638 


0.478 


0.762 


0.79 


0.36 


0.48 


0.30 


0.6 


0.537 


0.402 


0.746 


0.81 


0.57 


0.62 


0.33 


0.7 


0.431 


0.323 


0.721 


0.81 


0.63 


0.66 


0.30 


0.8 


0.320 


0.240 


0.688 


0.79 


0.65 


0.67 


0.23 


0.9 


0.200 


0.150 


0.643 


0.75 


0.64 


0.66 


0.15 


1.0 


0.070 


0.052 


0.588 


0.70 


0.61 


0.62 


0.05 



Table 1: The long-distance contributions to G2t(p), B^(p) and Bk{^) as determined 
using the ENJL model. Also shown are -Br-(i) using the one-loop short distance and B K ( 2 ), 
Br(2) usm S the two-loop short distance in Table 7. See Appendix A for the values of the 
parameters used. For the non-chiral cases one has to add 0.09±0.03 from the nonet vs 
octet difference, see text. 



(4.2) 

These integrals can be performed analytically but the result is rather cumbersome. 
The Euclidean continuation of q 2 we used is q\ = —q 2 . The result in the chiral limit 

becomes 

— 8B F 1 f 2 2 5 4 



{q 2 - m 2 K y 16tt 2 F 2 I ^ * 6 



- ~q* (4.3) 



and for q 2 = 



-8B 2 F 2 



(q 2 -m 2 K ) 2 16tt 2 F 2 



x { - 1(2^ - ml) („> - (2,4- - ml) h. 



2 ( 2 2 , I V + ™>k \\ 1 2 2 2i f ^ + m n 

mt 2 \ \ ml 



+m K [fi - m K In ( ^ ] ] - -m^ ( /i - In 




(4.4) 



These results allow to obtain the equivalent of (4.1) for the 0(p 4 ) coefficients. 



4.3 The Bk Parameter: Long Distance and Short Distance 

We now take the results from the ENJL evaluation of Ii K K °(q 2 ) both in the chiral 
limit and in the case of quark masses corresponding to the physical pion and kaon 
mass and use these to estimate B K and G 2 7- 

The final results for B K in the chiral limit, B^(fx) and G 2 7(/f) = AB^(fj,)/3 
are shown in Table 1. We have also shown the value of B K obtained there from 
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extrapolating the ENJL two-point function in the Euclidean domain to the kaon pole 
using Chiral Perturbation Theory, this is Bk(/j)- In the latter case we have to include 
the correction due to the difference between the octet and nonet case. This correction 
was estimated to be about 0.09 ± 0.03 in [2] and we take it as /i-independent. In the 
other columns in Table 1 various parts of the short distance correction are included. 
The realistic case, with non-zero quark masses in the long distance contribution 
to B K (fi), we have shown with the one- loop short-distance running, B K ^, two-loop 
short distance running with the scheme- dependence removed, Bf/^, as defined in eq. 
(A. 2), and the exact solution to the two-loop evolution equation with the scheme- 
dependence removed to the same order, B^F 2 ) as defined in Eq. (A. 3). For the latter 
short-distance contribution we have also shown the result in the chiral limit, B^ 2 y 
The rest of the parameters used are in App. A. Notice that the matching for all 
cases is acceptable. The quality of the matching for the real B K is as good as for 
d£{2) since they only differ by the /i-independent correction of 0.09 described above. 
So, in the chiral limit we get 

0.25 < B\ < 0.40 , (4.5) 
with non-zero quark masses we get 

0.50 < B^ onct < 0.70 (4.6) 

for the nonet case and 

0.59 <B K < 0.79 (4.7) 
for the real case. Notice that the large value of the chiral symmetry breaking ratio 

1.8 < |f < 2.4 (4.8) 

confirms the qualitative picture obtained in [2]. Finally, let us split the different 
contributions to the value of B K in the real case, 

B K = (0.33 ± 0.10) + (0.09 ± 0.02) + (0.18 ± 0.07) + (0.09 ± 0.03) (4.9) 

where the first terms is the chiral limit result, the second term are the 0(p A ) chiral 
logs at v — M p [2], the third term are the 0(p 4 ) counterterms and higher also at the 
same scale and the last term is the above mentioned contribution due to the rji — 1] 8 
mixing [2]. The error on the chiral log contribution is from varying v and the one on 
the counterterm contribution from looking at various different ways to extract the 
same counterterms as described in [2] plus some extra as an estimate of the model 
error. 

Notice that the last term in (4.9) is of order (m s — m d ) / 7V C 2 and it is not included 
in present lattice results. In fact, it introduces an unknown systematic uncertainty in 
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quenched and partially unquenched results which is difficult to pin down, see [11]. So 
the lattice results cannot be easily compared to ours. This term is also not included 
in determinations which are based in lowest order CHPT since it is higher order. 
Therefore again a direct comparison with our results has to be done carefully. 



5. AS = 1 Transitions: Long Distance 



In this section we use the AS = 1 two-point-functions H K ^ n ^(q 2 ) and H K ° n °(q 2 ) as 
defined in (3.1). We do not use the one with rjs since to order p 2 we do not get any 
more information out of that two-point function. It will provide extra information 
to 0{p A ) [12]. The result to lowest order in CHPT is given by 



n 



K + TT+ I „2 



(V 2 



AB 2 F 4 C 



IT^V 



(q 2 - m 2 K )(q 2 - m 2 ) 

2y / 2£ 2 F 4 C 
(q 2 - m 2 K )(q 2 - m 2 ) 



2 

-( 

3 



q 2 [Gs + -G 27 -2G',)+mlG' 



q 2 (-Gs + G 27 + 2G" 8 ) - m\G', 



(5.1) 



Here C of Eq. (2.7) has been chosen such that in the strict large- N c limit G 8 = 
G27 = 1. The coupling G' 8 is the coefficient of the weak mass term that does not 
contribute to K — > tttt at order p 2 but its value is important at 0(p A ) and higher 
and for some processes involving photons. The definition of the 0(p 2 ) Lagrangian, a 
discussion of the contributions from G' 8 and further references can be found in [12]. 

We have calculated the two-point functions in the chiral limit to extract the 
coefficient of q 2 and in the case of equal quark masses for an ENJL quark mass of 
0.5, 1, 5, 10, and 20 MeV in order to extract the coefficient of m 2 . 

As described in Section 3.2 we treat all coefficients (7* (//£,) as leading order in 
l/N c since they are enhanced in principle by large logarithms. We therefore obtain 
the matrix elements of Q 1 , Q 2 , Q3, Qa = Q 2 — Qi +Q3, Q5 and Q Q to next-to-leading 
order in 1/N C . 

5.1 Current x Current Operators 

The comments here are only valid for Qi, Q 2 , Q3, Q4, and Q$. The operator Q 6 is 
special and is treated separately in the next subsection. 

We can now use the method of [9] with the correct routing and obtain for the 
contributions to G 8 and G' 8 from Qi, Q 2 , Q3 and Q 5 [with = /! 2 /(167r 2 F 2 )]: 



G 27 (/i)[Qi] = GM\Q 2 ] = 1-3A, + 0(p 4 



GM\Q, 



l + ^A, + 0(p 4 ) 



9 



l + lA M + 0(^ 



G 8 (jm)[Q 2 ] - ■ 2 - /( 
GM[Q 3 ] = 2GM\Q 2 ] + 3G s (fi)[Qi] = + 0(p 4 ) , 
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G'M[Q 5 ] 



G'M[Qi] 



GM[Qs] 



+ O(p 4 ), 

o, 

^A, + 0( P 4 ), 

2G'MQ 2 ] = ^A, + 0(p i ), 
- 5 -A, + 0(p*). 



(5.2) 



Here and in the remainder G 8 {fi)[Qj\ stands for the long-distance contribution of 
operator Qi to G 8 when Ci(fx) is set equal to 1. The same definition applies to 
G' 8 (fx)[Qi] and G-n(n)[Qi\. In Tables 2 and 3 we dropped the argument (fx) for 
brevity. The results from the ENJL calculations are summarized in Table 2. The 
numbers in the columns 2 to 8 are always assuming [1 + (a s (fi)/n) r±j] Cj(fi) = 1 for 
the relevant operator. 

We get that G27(fi)[Qi] = G2i{^)[Q-2\ = G ? 27(a*) i n Table 1 and they are therefore 
not listed again. In addition all the other operators are octet so do not contribute 
to G27. We also have G' 8 (fi)[Qi] = 0, the operator Qi only contributes via _B^-like 
contributions which cannot have a contribution at q 2 = for equal quark masses 
since this type of contribution also produces G27 where such terms are forbidden. 
The approach to the chiral limit for the left-left current operators Qi, Q2, and Q 8 is 
such that the .B^Tike and Penguin-like contributions are separately chiral invariant. 
For the left-right current operator Q 5 this is not the case and it is only the sum of the 
.Bx-hke and Penguin-like contributions that vanishes for q 2 — > in the chiral limit. 
Notice that the results for small fx agree quite well with the results just using CHPT, 
eq. (5.2), but differ strongly for larger fx. The values (5.2) at fx — correspond to 
the factorizable contribution. 

We have also calculated the chiral logarithms that should be present in these con- 
tributions. Subtracting them made the extraction of the coefficient of m 2 to obtain 
G" 8 numerically much more convergent. The results for G 8 {fx)[Q^\ can be obtained 
from isospin relations from G$([/,)[Qi\ and G 8 (fi)[Q2\- The results for G 8 {fx)[Q§\ come 
from a large cancellation between the values of G 8 — 2G' % and G' 8 and have a somewhat 
larger uncertainty than the others. 

It should be noticed that in all cases the 1/N C corrections to the matrix elements 
are substantial. 

5.2 The Q 6 Operator: Factorization Problem and Results 

After Fierzing, the Qq operator defined in (3.5) 



gives both factorizable and non-factorizable contributions to the off-shell two-point 
functions to K° — > tt°, K° — > r/ 8 , and K + — > tt + . Here, we study for definiteness the 



Q 6 = "2E (1 + 75) q] (x) [q (1 - 75) d] (x) 



(5.3) 



16 



fi (GeV) 


G 8 [Qi\ 






s~i r /~\ "] 
G 8 [Q 5 \ 


GAQ-A 


G'slQa} 


G'siQs] 


0.0 


-0.667 


1.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.3 


-0.834 


1.271 


0.040 


-0.041 


0.070 


0.140 


-0.149 


0.4 


-0.930 


1.425 


0.060 


-0.109 


0.128 


0.256 


-0.297 


0.5 


-1.029 


1.600 


0.113 


-0.244 


0.206 


0.412 


-0.530 


0.6 


-1.130 


1.779 


0.168 


-0.460 


0.298 


0.596 


-0.868 


0.7 


-1.235 


1.962 


0.219 


-0.769 


0.399 


0.798 


-1.321 


0.8 


-1.347 


2.145 


0.249 


-1.178 


0.501 


1.002 


-1.908 


0.9 


-1.467 


2.325 


0.249 


-1.690 


0.598 


1.196 


-2.634 


1.0 


-1.597 


2.498 


0.205 


-2.308 


0.681 


1.362 


-3.504 



Table 2: The results for the long-distance contributions to Gs(^i) and Gg(/J,) from Qi to 
Qr> [Qi = Q2 — Qi + Q3] as calculated using the ENJL model via the two-point functions. 



K + — > 7r + two-point function, U K+7T+ (q 2 ) of Eq. (3.1). The factorizable contributions 
from Qq to this two-point function are 



- u p +K+ ( q )u P +n+ ( q )' . 

Here Cq^) is the Wilson coefficient of Q 6 , ITp(g) are two-point functions 



H%(q) = i J d 4 xe^ x (0\T(Pl(0)Pi(x))\0) = 



Z; 



q 2 — mf 



+ Z' 



(5.4) 



(5.5) 



with Pi(x) the pseudo-scalar sources defined in (3.2), and Ii p K- Sz2P ^+ {p^ q) the three- 
point function 



U P K s 3 ,P n+ip q) =fj d 4 x J d A ye i{ q .x- P .y) ( \T(P K -(x)S 32 (y)P 7T , 



(0))|0> 



with S 3 2(y) the scalar source 



S32 (y) = -[sd] (y) 



(5.6) 



(5.7) 



The last term in Eq. (5.4) corresponds to the diagram shown in Fig. 1(a). The first 
term is a contribution which is absent in the case of current x current operators. It 
is depicted in Fig. 5. 

In octet symmetry, to next-to-leading order we have [28] 

(0\dd + ss\0) 16 /n 9 9xr 4 „ . r rrx 3 5 . , 

9R J =l + 7^K+* + ^< (2L 8 + H 2 )--f^-3fj, K - -ix m (5.8) 
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Figure 5: The factorizable contribution for the Qq operator that is not well defined in the 
chiral limit. This contribution is not present for current x current operators. 



for the one-point function and 

n p *- 532 ^+(o, g ) = 



(q 2 - m 2 K )(q 2 - to 2 ) 



Bn 



X 



8 4 32 

1 + — 2 (2m 2 K + m 2 7r )(2L 6 - L 4 ) - (2q 2 + m\ + m\) L 5 + — 2 q 2 L 8 

*0 *0 *0 



1 q 2 



6 167T 2 F n 2 



In 



to 



K 



2 K - ml] 



m 2 In 



to 



m 



2 j+ m vs ln 

K J 



to; 



m 



K . 



mt 



+ 2^ + + 6^+12 lQn 2 F 2 



ln 



V 2 



m\ + ml 
167T 2 F n 2 8(ml 



K 



mt 



m„ , ln 



m 



m. 



— 3m w ln 



to; 



(5.9) 

*KJ \ m Kj\\ 

for the three-point function. Here and in the remainder the constants Li are defined 
at a scale u, Li = L\{u). and /ij = ln(TOj/V)/(167r 2 ) for i = n,K,i] 8 . 

At next-to-leading order, the expressions for the two-point functions were given 
for the octet symmetry case in [12]. So the second part in (5.4) can be written as 



n P K+K+ (q)U^ +7T+ (q) = 



yZx+Z^+ 2 R 2 

(q 2 - m 2 K )(q 2 - ml) 



x 



1 + (2to^ + to^)(4L 6 - L 4 ) + -^(ml + m 2 K )(4L 8 - L 5 ) 

4 7 5 5 

+-Ev( 2( l 2 ~ m l- m 2 K )(2L 8 - H 2 ) - -fi* - -fi K - 



F 2 



12' 



(5.10) 



Therefore this order, in octet symmetry, the factorizable contributions to H K+W+ (q) 



from Qq are 



n 



K + TT + 

QQ,Fact 



V z K +z n + 



{q 2 - m 2 K )(q 2 - ml) 



C AS =iC 6 (fx) 16B 2 (fi) 



x 



2q 2 {L 5 - {2L 8 + H 2 )} + m\{2L 8 + H 2 ) - —fi K + — a m - 



12' 



16' 



16' 



18 



24 16tt 2 F 2 



In 



m 



K 



2{m\ 



ml 



In 



m 



+ 



ml) 



K 



48 16tt 2 F 2 \v 2 ) 
m\ + m 2 1 
16tt 2 F 2 32{m 2 K -m\) 



m 2 In 



T?8 



to; 



to 



A" , 



m 



K , 



3m n ln 



to: 



m 



K , 



m 



V8 



m 



K , 



(5.11) 



As it is well known the order p° contribution from Qq vanishes [29] and the first 
non-trivial contribution from this operator is of order 3 p 2 . This happens here as an 
exact cancellation between the two types of factorizable contributions at order p°. 
As a result there is a very large cancellation between the two types of factorizable 
contributions at order p 2 . We get 



Q 6 ,Fact 



16 CM 



BM 



U- 



3 1 

T6167T 2 



(5.12) 



and 



G' 



Q 6 ,Fact 



8 cm 



B 2 M 
F 2 



(2L S + H 2 )-^- 



24 16tt 2 



21n 



(v) 



+ 1 



(5.13) 



The mass m L above has to be understood as an infrared cut-off as we have done 
the chiral limit m^ = m^ = mx = rn m — > 0. The factorizable contribution to G$ 
and Gg from Qq is therefore not well defined. It has an infrared divergence. The 
divergence is related to the divergence in the pion scalar radius in the chiral limit. 
Since Qq is an 8^ x l fi operator we know from CHPT in the non-leptonic sector 
that to lowest order in the counting there, no infrared divergences are present in the 
two-point function U K+7T+ (q 2 ). These infrared divergences are therefore spurious and 
must be cancelled by another contribution. The only possibility is that it cancels out 
with the non-factorizable contribution also coming from Qq. We will see below that 
this is indeed the case. Notice also that since G 8 and G' 8 are 0(p 2 ) couplings, Eqs. 
(5.12) and (5.13) are exact for the factorizable contributions. 

Unfortunately, the non-factorizable contributions can only be calculated at present 
in a model dependent way. In the 1/N C expansion, the infrared divergent part of G§ 
and G' s can in fact be calculated analytically using the 0(p 2 ) CHPT Lagrangian. We 
can therefore subtract it. It follows from the diagrams shown in Fig. 3, (b),(c), (e), 
and (f) and by using CHPT for the A-boson vertices which is valid for small ji. For 



equal masses m 2 K = m 2 



m 



m 



2 T we obtain 



n 



2B 2 F 2 



(q 2 -m 2 K )(q 2 - m 2 ) 



C AS=1 C 6 (fi)4B 2 M 



3 The order p 2 chiral logs were called order p°/N c contributions in [30]. 
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(2tt)< 



({r e + Qe) 2 + m 2 L ) (r| + m|) 



The non-factorizable (NF) part above in the limit — * leads to 



G, 



Q 6 ,NF-0(p 2 ) 



16 C 6 (/i) 



B 2 M 3 1 



Fo 2 



16 16tt 2 



2 In 




(5.14) 



(5.15) 



and 



G' 



Q 6 ,NF-0(p 2 ) 



8 cm 



B 2 M 5 1 
F 2 24 16tt 2 



2 In 



+ 1 



(5.16) 



There is a very large cancellation between the factorizable parts in (5.12) and 
(5.13) and the non-factorizable part in (5.15) and (5.16) both for the IR divergent 
part and for the large 1/N C constant part. Summing up the exact factorizable result 
and the infrared divergent non-factorizable part we get 



G, 



Qe, 0{p 2 ) 



and 



G' 



Qe, 0{p 2 ) 



16 CM 



B 2 M 
Fo 2 



LM - 



IGtt 2 VS '"(!!.) " 96 



(5.17) 



8 CM 



B 2 M 
Fi 



[2L 8 + H 2 )(u)- 



12 16tt 2 



In 



// 



(5.18) 



It is then a non-trivial check of the validity of the model used that the non-factorizable 
part indeed contains the correct infrared logarithms needed to cancel the factorizable 
ones. The ENJL model used here does. 

Notice in (5.17) and (5.18) all the dependence on the IR scale , m|, drops out 
as it should and the scale in the logarithm becomes \n(n/v). So in the chiral limit 
and next-to- leading in 1/N C , the scale dependence on the short-distance scale gets 
compared to the scale where the CHPT constants are defined. 

The result above shows that at least the B 6 parameter defined as usual as the 
ratio of the non-factorizable contributions over the vacuum saturation result (VSA) 
is not well defined. It is therefore necessary to give another definition for this B 
parameter. The cancellation of the infrared divergence found here is probably also the 
source for the large cancellations found between the factorizable and non-factorizable 
contributions in earlier work. Notice also that the 1/N C finite term in (5.12) is larger 
than the leading in 1/N C result and with opposite sign. It is clear that it can be 
dangerous not to have and analytical cancellation of both the IR divergent part and 
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the 1/N C constant as we have. This can explain also some discrepancies for the B 6 
parameter results in the literature, B$ is just not well defined. 

The way we treat our results is that we remove the exact infrared logarithm 
from our ENJL calculation by adding equations (5.12) and (5.13) which are exact 
and model independent to the ENJL results. In this way we also remove the IR 
divergence of the non-factorizable part exactly. We chose the reference scale \i = M p 
to do the subtraction. We generate the mass m 2 L by putting small current quark 
masses. The remaining factorizable factor, i.e. the part from the constants L 5 , L 8 , 
and H 2 are then evaluated at a scale v = M p . This corresponds for the leading in 
1/N C contribution to G$ and G' 8 from Qq 

= (-9 ± 14) CM (5.19) 
using 

L 5 (M p ) = (1.4 ±0.3) • 10~ 3 
(2L 8 + H 2 )(M p ) = (0.7±1.1)-10- 3 . (5.20) 

We have used here the value of B and F from the ENJL model. The value of 
2L 8 + H 2 is derived from the canonical value for L 8 (M P ) = (0.9 ± 0.3) • 10~ 3 and the 
value for (2L 8 - H 2 )(M p ) = (2.9 ± 1.0) • 10~ 3 from [33]. The large error for G' S ENJL 
in (5.19) is because of the large cancellation in the value for 2L 8 + H 2 . Notice that 
the size of the subtracted terms in Gf NJL is about +40C 6 (/x) for m 2 L = m^rriK and 
varies very fast with uil- 

Our calculation agrees with the one of [30] when the appropriate identifications 
are made. The large cancellation between the factorizable and non-factorizable parts 
where also observed there. They were however not identified as an exact cancellation 
of infrared divergences. In fact, at the order the calculation was done in [30] the 
cancellation of the l/N c factorizable and non-factorizable pieces is very large, and in 
their language 4 one should get B 6 very near to one. They get indeed B e very close 
to one. 

The non-factorizable non-divergent part has corrections from higher order terms 
in the chiral Lagrangian which we calculate numerically using the ENJL model. We 
have included them and these give therefore the numerical differences between our 
results and the ones in [30]. 

Before we present the results for G 8 (/i) and G' 8 (fx) from Q 6 from our ENJL 
calculation we need to include one additional remark. The vector and axial-vector 
currents used in the previous section are uniquely identified both in the ENJL model 
and in QCD. There is however no guarantee as remarked in [23] that the same is 
true for the scalar and pseudo-scalar densities. Here we renormalize the ENJL scalar 

4 As we said B 6 is not well denned. We come back to this question in Section 7. 



G 



ENJL 



7V„ 



[-38 ±8) CM and G 



1 ENJL 
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fj, (GeV) 










GslQe] 


G'siQe] 




ENJL 


ENJL 


(i) 


(i) 


(2) 


(2) 


0.3 


-118 


-69 










0.4 


-103 


-53 










0.5 


-93 


-41 


-21.1 


-9.3 


-6.4 


-2.8 


0.6 


-88 


-32 


-23.9 


-8.7 


-14.7 


-5.3 


0.7 


-84 


-25 


-25.9 


-7.7 


-20.1 


-6.0 


0.8 


-82 


-20 


-27.9 


-6.8 


-24.5 


-6.0 


0.9 


-82 


-17 


-30.0 


-6.2 


-28.4 


-5.9 


1.0 


-83 


-15 


-32.4 


-5.9 


-32.4 


-5.9 



Table 3: Results for the long-distance contributions to Gs and from Qq as calculated 
using the ENJL model via the two-point functions for the non-factorizable part and adding 
the model independent factorizable part in (5.12) and (5.13). The last 4 columns include 
the renormalization of scalar and pseudo-scalar densities to one-loop W and two-loops ^ 
in QCD. The short-distance anomalous dimensions for Bo(/j,) at scales below 0.5 GeV blows 
up. 



S(x) and pseudo-scalar P(x) densities by the values of the quark condensates in the 
chiral limit: 

There is an analogous equation for the pseudo-scalar density. This factor should be 
remembered when using the Wilson coefficients from our results. The values we have 
used are B$ CD (lGeV) = (1.75±0.40) GeV in the WS scheme [33, 34], and B$ NJL = 
2.80 GeV [21]. We have also included the QCD scale dependence of the B parameter 
to two-loops. We show in Table 3 the results for G 8 (/•<.) [£*?6] an< ^ G' 8 ([i)[Q e ] without 
the renormalization factor of Eq. (5.21), columns labelled ENJL, and including the 
renormalization factor of Eq. (5.21) both to one-loop, columns labelled ^\ and two- 
loops in QCD , columns labelled ( - 2 \ Notice B (fi) = —(qq)(fi)/FQ and this factor is 
responsible for most of the running of (5e[31]. 



6. The Order p 2 Full AS = 1 Couplings 

We use here the results of [7] and [6] for the AS = 1 QCD anomalous dimensions to 
one- and two-loops respectively to obtain final values. The solution for the Wilson 
coefficients are given in [7, 19] at two- loops using an expansion in a s . Whenever the 
values of Aqcd are needed in the MS scheme with three flavours we use the expanded 
in a s formulae [5] from a s (M T ) = 0.334 ± 0.006 with M T = 1.77705 ± 0.00030 GeV 
[5] and get A^ D = 220 MeV to one-loop and A ( q ] cd = 400 MeV to two-loops. The 
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H (GeV) 


G27 


G$ 


G's 


0.5 


0.399 


4.45 (4.55) 


0.739 (0.761) 


0.6 


0.351 


4.26 (4.34) 


0.686 (0.710) 


0.7 


0.291 


4.21 (4.28) 


0.703 (0.727) 


0.8 


0.221 


4.25 (4.30) 


0.767 (0.789) 


0.9 


0.141 


4.33 (4.37) 


0.847 (0.866) 


1.0 


0.050 


4.44 (4.46) 


0.923 (0.935) 



Table 4: The final results for the three 0(p 2 ) couplings using the one- loop Wilson coeffi- 
cients. The numbers in brackets refer to using Q±, Q 2 , and Qq only. 

values of the Wilson coefficients we use for AS = 1 [7, 19] and for AS = 2 [32] are 
in the Appendix. We also include there the scheme dependent constants r\ needed 
for the two-loops short- distance running in the NDR scheme we use. 

We now show in Tables 4 and 5 the results for the coefficients G 27 , G 8 and G' s . 
The numbers in brackets refer to keeping only Qi, Q2, and Q 6 . Most of the difference 
is due to Q4. 

The matching for the one-loop running of the Wilson coefficients is very good. 
We obtain a value of G$ ~ 4.3 and G' 8 ss 0.8. If we look inside the numbers, for 
G$ the contribution via Q\ is fairly constant over the whole range but there is a 
distinct shift from Q 2 to Qe for lower values of /i. The operator Q 2 remains the 
most important over the entire range of fi considered. For G' s similar comments 
apply except that Q\ doesn't contribute. Typically G27 is somewhat low compared 
to the experimental number and we have not as good matching as in the octet sector. 
Notice though that it gets somewhat more stable in the range between 0.5 and 0.8 
GeV as one expects from the validity of the low-energy model. 

When two-loop running is taken into account in the NDR scheme the numbers 
do not change so much. The effect of the r\ constants in this scheme is however very 
large and causes a significant shift in the numbers. 

The numbers for the octet case are somewhat stable in the range \i = 0.8 to 1.0 
GeV but there is where the ENJL model is expected to start deviating from the true 
behaviour. 

Notice that at large N c , G$ and G27 are both 1. Adding 1/N C corrections G27 
decreases by a non-negligible factor around two to three, while the G$ coupling gets 
enhanced up to G§ = 6.2 ± 0.7. The short-distance enhancement is almost a factor 
of two for the whole range of fi. The rest of the enhancement, namely a factor two 
to three is mainly due to the large value of the long-distance contribution to the 
Penguin-like coupling c. The bulk of the long distance part enhancement of the 
coupling c comes from Q 2 and Qq. There is also a small contribution to G$ in the 
right direction from the S^-like coupling b from both Q 2 and Q\. 
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H (GeV) 


G 2 -j 


G% 




0.5 


0.182 


11.20 (12.4) 


1.60 (1.75) 


0.6 


0.249 


7.30 (7.8) 


1.13 (1.22) 


0.7 


0.230 


6.30 (6.6) 


0.99 (1.10) 


0.8 


0.184 


5.88 (6.2) 


0.97 (1.08) 


0.9 


0.121 


5.73 (5.9) 


0.99 (1.11) 


1.0 


0.044 


5.61 (5.8) 


1.03 (1.14) 



Table 5: The final results for the three 0(p 2 ) couplings using the two-loop Wilson coeffi- 
cients with the inclusion of the r\ factors. The numbers in brackets refer to using Q±, Q 2 , 
and Qq only. 



H (GeV) 


One-Loop 


Two-Loops 


0.5 


14.3 


78.5 


0.6 


15.6 


37.5 


0.7 


18.6 


35.0 


0.8 


24.6 


40.8 


0.9 


39.2 


60.1 


1.0 


113.2 


162.4 



Table 6: The final results for the ratio \A /A 2 \ to 0(p 2 ) using the one-loop short-distance 
running and the full scheme independent two-loops short-distance running. 

The final results for the ratio \A /A 2 \ at 0(p 2 ) (2.10) are in Table 6. The stability 
we get for the one-loop short- distance is not bad, and there is some minimum around 
0.7 GeV for the two-loop running. We get in general too large values for this ratio 
compared to the experimental 16.4 value (2.5) due to the somewhat small value of 
G27 we get. 

In order to show the improvement with previous results and the quality of the 
matching we have shown in Figure 6 for G 27 (fx) the lowest order result Eq. (5.2), the 
ENJL result for the same quantity and the final result for G 27 with the two-loop short 
distance included. We have similarly plotted Gs[Qi](fJ>) and Gs[Q 2 ](fi) both from the 
lowest order result Eq. (5.2) and from the ENJL model. We also showed the full result 
for G' 8 when the two-loop running is included properly. Similar improvements of Eq. 
(5.2) and (5.17) can be seen by plotting the other results with the corresponding 
ones from Tables 2, 3, and 5. 

In summary, the results we get for G$, G 2 j, and G' s are 

4.3 < G 8 < 7.5 
0.8 < G' 8 < 1.1 
0.25 < G 27 < 0.40 (6.1) 
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Figure 6: The improvement of the behaviour with \i of several quantities. Shown are the 
lowest order result, ENJL result, and when short-distance running is added. 



The bounds have been chosen by looking at both the one-loop and two-loop results 
in the stability regions in Tables 4 and 5. From (6.1) we can extract the values 



-0.75 < b < -0.50 
1.7 < c < 3.7 

and we have fixed a = 1 as explained before. For the AI = 1/2 rule we get 

\A W 



15 < 



A? 



< 40 



(6.2) 



(6.3) 



to order p 2 . 

We get a huge enhancement due to the c-coupling, it is therefore interesting 
what do other calculations predict for this coupling. One model where this coupling 
can be easily extracted is the effective action approach [35]. To order 1/N C one gets 
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[35, 36] 



c = CM - 1 + KeCtiti) + C 2 {ll) + L 10 )(/i) 



UeCM + CM (2H, + L 10 )(/i) 



b 2 M 

+ 0(l/iV c 2 ) (6.4) 



with li = M p , a s {M p ) = 0.70, B (M P ) = 1.4 GeV, and (2#i + L 10 )(M p ) = -0.015 
[21], we get 

c = 0.95 ±0.40. (6.5) 

The reason why c is smaller than the present work results is that the long-distance 
mixing between Q 2 and Qq is not well treated in this model. In fact this contribution 
is model dependent already at 0(1/N C ). For instance, it appears in terms of the 
short-distance value a s (M p ). It is clear that at such scales one has to treat the 
long distance contributions in a hadronic model and the a s (M p ) above will appear 
enhanced. Nevertheless, the extra contribution to c coming from the operator Q 2 
[35, 36] both from short- distance origin, namely the term C 2 (/jl) — 1, and from long- 
distance origin, namely the part proportional to 2Hi + L w , give some insight on the 
potentially large value of c. 

We cannot easily compare our result with those of [37], their method of calcu- 
lating the low energy part has no obvious connection to the short-distance evolution 
and their results cannot be directly compared to ours. The results from the lattice 
[11] are at rather high values of the quark masses and can thus also not be simply 
compared to our results. 

As stated above, we agree with the calculations of [30] for low values of the 
scale ix where we should agree but deviate significantly at higher scales. The earlier 
Dortmund group results [38] are thus also expected to have significant corrections. 
The attempts at calculating via more inclusive modes [39] have very large QCD 
corrections [35, 40]. We see the remnant of this in the large corrections from the 
T\ terms, see Appendix A. The short-distance factors are in fact one of the bigger 
remaining sources of uncertainty. 



7. Results and Conclusions 

The main results of this paper are the results for the 0(p 2 ) couplings G 8 , G 27 , and 
G' 8 as a function of cut-off ix for the various operators Qj, j = 1, • • • , 6, as given in 
Tables 1, 2, and 3. In addition we have corrected our earlier results for B k (li) for 
the routing problem as described in Section 4.1 and presented those in Table 1 as 
well. 
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The other main result of this paper is the observation that in the chiral limit the 
factorizable contribution from Qq is not well defined due to an infrared divergence 
and we expect that similar problems will show up for the current-current operators 
when we try to calculate higher order coefficients in the weak chiral perturbation 
theory Lagrangian. We showed that the total contribution of Qq obtained after 
adding the non-factorizable and factorizable parts is however well defined. We also 
expect that he same solution will hold for coefficients of higher order operators in 
the chiral lagrangian. A corollary of this observation is that the use of 5-factors in 
the chiral limit as is common in other treatments of weak non-leptonic operators is 
not possible in the way they are defined, namely the whole result normalized to the 
VSA result. 

One could use the leading result in l/N c as an appropriate starting point for 
normalizing the 5 6 -parameter in the chiral limit, this is keeping only the L 5 , L 8 , and 
H 2 terms but this is difficult to implement for lattice gauge theory calculations. In 
fact, what in practice people have used [19, 30] for the VSA, i.e. the factorizable 
part of Qq, has been just the large N c part. Of course, this is not in agreement with 
what is done with other B parameters for current x current operators like B^ where 
the 1/N C factorizable part is always included in the VSA result. After the problems 
we encountered and the importance of the B parameters to normalize results from 
different techniques, we believe a new consistent definition of the B parameters should 
be looked for or just abandon the use of B parameters and quote matrix elements 
values. We also emphasize that caution should be taken when combining results from 
different methods for the factorizable and non-factorizable contributions. 

When we combine our main results with the Wilson coefficients at one-loop we 
get nice stable results. Using the Wilson coefficients at two-loops with the inclusion 
of the r\ factors which as we argued in Section 3.2 is necessary we obtain relatively 
stable values for G$ and the coefficient of the weak mass term G' 8 with 

4.3 < G 8 < 7.5 
0.8 < G' 8 < 1.1 
0.25 < G 27 < 0.40 (7.1) 

The main uncertainty here is in fact coming from the short-distance coefficients for 
the octet case and from the long-distance for the 27-plet case. For the G27 coupling 
we obtain a somewhat small value compared to the experimental one. This translates 
into the following results for the A/ = 1/2 rule in the chiral limit 

(2) 

<40. (7.2) 

These results are somewhat large. Nevertheless, we would like to emphasize that we 
have obtained these results from a next-to-leading in 1/N C long-distance calculation 



15 < 



A, 
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and we have passed from the large N c result | ^4 / ^4 2 1 tv c — V% to values around 20 to 
35. One can certainly expect non- negligible 1/N% corrections to our results but the 
huge enhancement is there. We would also like to stress that we have no free input 
in our calculation. All parameters have been determined from elsewhere. 

From the results above we have also obtained the couplings b = G27 — 1 and 
c = (3G 8 + 2G 27 )/5 - 1 

-0.75 < b < -0.50 

1.7<c<3.7. (7.3) 

Here is then one of our main results, the AI = 1/2 rule enhancement comes from 
the Penguin-like topologies (c) in Figure 1, both from Q 2 which dominates for high 
values of ji and from Q 6 which dominates for small values of /i. 

In addition we obtain a value for chiral limit value B\ as defined in Eq. (A. 2) 

0.25 < B\ < 0.40 (7.4) 

and the value for the Bk parameter in the real case 

0.59 < B K < 0.79 . (7.5) 

These two results confirm the ones in [2]. Notice that the different short-distance 
contribution from M w until the charm quark mass to G27 and B\ has produced 

tt- * 11 (7-6) 

instead of 3/4. 

So we have obtained quite good matching for G%, G' 8 , and B K for values of /i 
around 0.7 — 1.04 GeV and for G27 for values of fi around 0.6 — 0.8 GeV. We obtained 
values for the three parameters of order p 2 not too far from the experimental ones 
and a quantitative understanding of the origin of the AI = 1/2 enhancement. Notice 
that the values of the cut-off we use to predict our results are not extremely low as 
in other 1/N C approaches, still one would like the matching region to be larger and 
for somewhat larger values of the cut-off. 
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A. AS = 1 and AS = 2 Wilson Coefficients 



In this section we give the numerical values of the Wilson coefficients for the basis 
of AS = 1 operators in (3.3), (3.4), (3.5), and for the AS = 2 operator in (3.7). We 
give them for the relevant values of the renormalization scale. We have extensively 
used the formulae in [19]. 

In all cases we have used a s (M\y) = 0.121 ± 0.002, obtained from LEP measure- 
ments at the Z-peak [5] and then run to two loops to M w = (80.41 ± 0.10) GeV, 
ct s {^b) = 0.232±0.003 obtained from QCD sum rules in the T system [41], a s (M T ) = 

0. 334 ± 0.006 [42], we then run this value to two-loops up to m c (m c ) = (1.23 ± 0.05) 
GeV in the MS scheme [43] and get a s (m c ) = 0.42 ± 0.01. In our approach [19], 
the penguin operators only get generated from the charm quark mass down since the 
very small part due to the top quark is not relevant here. We have used the exact 
solutions of the renormalization group for the running of a s and the quark masses. 
For C3, C4, C5, and Cq we set the small imaginary part due to the top loop to zero. 
The results to one-loop accuracy are in Table 8 and for two-loops are in Table 9. In 
the two-loop case we are using the NDR scheme. 

If we run the short-distance contribution to two-loops, the matching condition in 
(3.14) sets a further coefficient which rends the matrix element scheme independent, 

1. e. 

l + ^n. (A.l) 

7T 

In the case of the AS = 2 operator, the full two-loop calculation can be found in 
[32]. In the NDR scheme we have 77 = —7/6 for the AS = 2 operator. We obtained 
it from the right eigenvalues of r\ in [19]. So we define 



Bk = (l + a " iP) 



7T 



. 72 P2I1 



M/i)p //31 BM . (A.2) 



With (3i = -9/2, fo = -8, 71 = 1, and 72 = -17/48. From the discussion in [19] 
and Section 3.2 it can be seen that this definition is scheme and renormalization 
scale independent. We have shown in Table 7 this factor in front of B K (fj.) for the 
case of one-loop running labelled One-loop, two-loop running with r\ = 0, labelled 
Two-loops, r\ at its value, labelled Scheme-Independent (SI), and a version where we 
use the exact solution of the two-loop running with 77 included so as to cancel the 
full scheme-dependence there too, i.e. 



S*=[1 + !^J [a s (^B M . (A.3) 



This is labelled exp in Table 7. 

The short-distance results for the AS = 1 Wilson coefficients to two-loops and 
including the (A.l) term which can be found in the NDR scheme in [7, 19] for 
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MGeV) 


One-Loop 


Two-Loops 


SI 


exp 


0.50 


1.04 


1.11606 


0.46894 


0.63252 


0.60 


1.08 


1.14653 


0.76287 


0.82690 


0.70 


1.12 


1.18208 


0.87769 


0.91869 


0.80 


1.15 


1.21045 


0.94751 


0.97817 


0.90 


1.17 


1.23348 


0.99680 


1.02160 


1.00 


1.19 


1.25267 


1.03445 


1.05546 



Table 7: The coefficients to transform Bk(/j) into Bk- See text for an explanation of the 
different columns. 



/i(GeV) 


Ci 


c 2 


c 3 


c 4 


c 5 


c 6 


0.50 


-0.96466 


1.59028 


0.01647 


-0.03796 


0.01116 


-0.04663 


0.60 


-0.84146 


1.49560 


0.01067 


-0.02626 


0.00801 


-0.03037 


0.70 


-0.75899 


1.43423 


0.00710 


-0.01839 


0.00576 


-0.02039 


0.80 


-0.69875 


1.39058 


0.00468 


-0.01263 


0.00403 


-0.01356 


0.90 


-0.65222 


1.35759 


0.00292 


-0.00816 


0.00264 


-0.00854 


1.00 


-0.61482 


1.33159 


0.00158 


-0.00455 


0.00149 


-0.00467 


Table 8: Wilson Coefficients of the operators Q\ 


to Qq at one- loop. 


/i(GeV) 


Ci 


c 2 


c 3 


c 4 


c 5 


c 6 


0.50 


-0.80875 


1.48719 


0.13750 


-0.26345 


0.01338 


-0.27035 


0.60 


-0.74066 


1.43763 


0.05198 


-0.11330 


0.02483 


-0.09696 


0.70 


-0.65083 


1.36940 


0.03088 


-0.07225 


0.02160 


-0.05673 


0.80 


-0.58661 


1.32243 


0.02097 


-0.05124 


0.01849 


-0.03770 


0.90 


-0.53854 


1.28836 


0.01516 


-0.03796 


0.01595 


-0.02631 


1.00 


-0.50087 


1.26236 


0.01133 


-0.02861 


0.01388 


-0.01860 



Table 9: Wilson Coefficients of the operators Q\ to Qq at two-loops in the NDR scheme. 

instance, are in Table 10. Here we give the one-loop results in Table 8, two-loop 
results with 5 r\ — at two-loops in Table 9 and the one with the scheme dependence 
properly removed, including r±, in Table 10. It can be seen that the change from one 
to two-loops in the NDR scheme is not so large but inclusion of the r\ makes a large 
change. 



5 Of course all quantities here are matrices. 



30 



MGeV) 


Ci 


c 2 


c 3 


c 4 


c 5 


c 6 


0.50 


-3.73959 


4.02465 


0.31282 


-0.43205 


0.03267 


-0.33360 


0.60 


-1.89282 


2.35657 


0.10089 


-0.16996 


0.02789 


-0.10140 


0.70 


-1.41708 


1.95062 


0.05666 


-0.10722 


0.02258 


-0.05469 


0.80 


-1.17990 


1.75588 


0.03741 


-0.07706 


0.01881 


-0.03390 


0.90 


-1.03270 


1.63865 


0.02665 


-0.05877 


0.01597 


-0.02190 


1.00 


-0.93034 


1.55917 


0.01979 


-0.04625 


0.01374 


-0.01397 



Table 10: Wilson Coefficients of the operators Q\ to Qq at two-loops, the NDR scheme 
dependence is removed as in discussed in Section 3.2. 
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